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^ ! ABSTRACT: We study the time evolution of the reduced density operator for a class 

K^ ' of quantum Brownian motion models consisting of a particle moving in a potential V{x) 

^ . and coupled to an environment of harmonic oscillators in a thermal state. Our principle 

Q I tool is the Wigner function of the reduced density operator, and for linear systems we 

[^^ . derive an explicit expression for the Wigner function propagator. We use it to derive 

O i two generalized uncertainty relations. The first consists of a sharp lower bound on the 

^ ' uncertainty function, U = (Ap) (Ag) , after evolution for time t in the presence of an 

->«^ ! environment. The second, a stronger and simpler result, consists of a lower bound at time 

qh| t on the quantity, A = U — Cp^, where Cpq = ^{ApAq + AqAp) . {A is essentially the area 

^ I enclosed by the 1 — cr contour of the Wigner function) . In both cases the minimizing initial 

bS)- state is a correlated coherent state (a non-minimal Gaussian pure state), and in the first 

j>! I case the lower bound is only an envelope. These generalized uncertainty relations supply 

•^ ■ a measure of the comparative size of quantum and thermal fiuctuations. We prove two 

r> I simple inequalites, relating uncertainty to von Neumann entropy, — Tr(plnp), and the von 

c^ I Neumann entropy to linear entropy, 1 — Trp . We also prove some results on the long-time 

limit of the Wigner function for arbitrary initial states. For the harmonic oscillator the 
Wigner function for all initial states becomes a Gaussian at large times (often, but not 
always, a thermal state) . We derive the explicit forms of the long-time limit for the free 
particle (which does not in general go to a Gaussian) , and also for more general potentials 
in the approximation of high temperature. We discuss connections with previous work by 
Hu and Zhang and by Paz and Zurek. 
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l.INTRODUCTION 

Quantum Brownian motion models have been the subject of a number of studies of 
many years [1,2,3,4,5,6,7,8,9,10,11]. Many reasons for the interest in these models may be 
found: they permit the possibility of studying in some detail the approach to equilibirum 
in non-equilibrium system; they arise in studies of macroscopic quantum effects; and they 
are related to the question of dissipation in tunneling. Most recently, they have been 
studied in the contexts of quantum measurement theory, decoherence, and the quantum 
to classical transition. It is these contexts towards which the present work is directed. 
We are, in particular, interested in the emergence of classical behaviour in open quantum 
systems. 

The quantum Brownian motion models belong to an important class of non-equilibrium 
systems in which there is a natural separation into a distinguished subsystem, S, and the 
rest (the environment). The distinguished subsystem S is often referred to as an open 
quantum system. One is interested in the behaviour of 5, but not in the detailed behaviour 
of the environment. S is most completely describe by the reduced density operator p 
obtained by tracing out over the environment states. One's goal is then to obtain an 
evolution equation for p, which will in general be non-unitary, from which one may calculate 
the probabilities of any observables referring to S only. In the one-dimensional QBM 
models studied in this paper, S consists of a particle of mass M moving in a potential V{x) . 
S is linearly coupled to an environment, consisting of a large number of harmonic oscillators 
in a thermal state at temperature T, and characterized by a dissipation coefficient 7. 

Given such a model, there are many interesting questions one can then ask about it. 
Under what conditions is there suppression of interference between localized wave packets? 
Under what conditions does the Brownian particle evolve approximately classically? How 
big are the fluctuations about classical predictability? Are these fluctuations larger than 
the inescapable quantum fluctuations? Is there a generalization of the uncertainty principle 
to include environmentally-induced fluctuations, representing the smallest amount of noise 
the system must suffer? What sort of states are most stable in the face of these fluctuations? 
Does entropy /uncertainty increase as time evolves? Towards what sort of states does the 
system evolve in the long-time limit? Does the system tend towards thermal equilibrium 
in the long-time limit? 

2 



Many of these questions have been addressed before [12,10,13,14,15,11,16,17] and some 
of them wiU be the topic of the present paper. We will address in particular, the question 
of the generalization of the uncertainty principle, and the question of the long-time limit 
of arbitrary initial states. The methods we will develop, however, will be applicable to 
most of the other questions listed above. 

Our starting point is the observation that all of the above questions depend on the 
dynamical evolution of the system. It is therefore of crucial importance to possess a clear 
picture of that evolution. This we obtain by focusing on the evolution of the Wigner 
function of the reduced density operator. In particular, the tool we found to be of greatest 
use is the Wigner function propagator. The reason why the Wigner representation is so 
clarifying is that, at least for linear systems, the unitary part of the evolution corresponds to 
transporting the Wigner function along the classical phase space trajectories. The Wigner 
function propagator in this case is then just a product of delta-functions. The same holds 
to leading order in an h expansion for non-linear systems. It is therefore possible to cleanly 
separate the unitary effects from the non-unitary effects induced by the environment. 

In the presence of an environment, the Wigner function propagator may be calculated 
exactly for linear systems, and it is just a Gaussian. It permits the explicit calculation of 
all moments of p and q at any time for arbitrary initial states in terms of the moments 
at the initial time. For systems with more general potentials, we cannot calculate it 
exactly, but some interesting aspects of the evoluion may be extracted using a semiclassical 
approximation. 

We describe the features of quantum Brownian motion models in Section II, and we 
calculate the Wigner function propagator for linear systems in Section III. 

In Section IV, we use the results of Section III to derive generalized uncertainty re- 
lations for quantum Brownian motion models. Two relations are derived. The first is a 
lower bound over all initial states on the uncertainty function 

U = {ApfiAqf (1.1) 

at an arbitrary time t. This lower bound represents the least amount of noise, both 
quantum and environmentally-induced, the system must suffer after evolution for time t in 
the presence of an environment. The lower bound is not the time evolution of a particular 
initial state, but is actually an envelope - the initial state achieving the lower bound at 
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time t is different for each time. The second relation is a lower bound on the related 
quantity, 

A^ = (ApfiAqf - hApAq + AqApf (1.2) 

where Aq = q ~ (q) and likewise for Ap. A is essentially the area enclosed by the 1 — 
a contour of the Wigner function. A slightly stronger version of the usual uncertainty 
principle exists in terms of A, taking the form A > h/2 [18]. The time evolution of A 
turns out to be much simpler than that of U, and the lower bound on it at any time is 
the time evolution of a particular initial state, hence this result is stronger than the first 
one. In both case the minimizing initial state is a correlated coherent state. A non-minimal 
Gaussian pure state, of the form ijj(x) = e^^"'^^ '^ , where a, b are real. (This is clearly 
some kind of squeezed coherent state, but we choose to follow the nomenclature of Ref.[18]). 
The detailed forms of the lower bounds on U and A are discussed in Section V, in the case 
of an ohmic environment. 

In Section VI we discuss the connection of these measures of uncertainty with von 
Neumann entropy, S[p] = — Tr(plnp). We prove an inequality relating uncertainty to von 
Neumann entropy, which in the regime of large uncertainty has the form, U > A > he. 
We also exhibit the connection with the linear entropy, 5"^ = 1 — Trp . 

In Section VII, we consider the long-time limit of an arbitrary initial state. For the 
harmonic oscillator in an ohmic environment, all initial states go to a Gaussian Wigner 
function in the long-time limit, which is not always a thermal state. In Section VIII we 
show how some of our results may be generalized to systems with more general potentials. 
We also consider the case of the free particle. 

We summarize and discuss our results in Section IX. 



2. QUANTUM BROWNIAN MOTION MODELS 

We are concerned in this paper with the class of quantum Brownian models consisting 
of a particle of large mass M moving in a potential V{x) and linearly coupled to a bath 
of harmonic oscillators. (This section is entirely review of standard material [19,4,8,20].) 
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The total system is therefore described by the action, 



STot[x{t),qn{t)] 



dt 



1 



™2 



Mx^ - V{x) 



+Y.J dt 



-rUnql - -rUnUjlql - CnQnX 



(2.1) 



Quantum-mechanicaUy, the total system will be described most completely by the density 
matrix, pt{x , qn, x' ^ q'^ . However, we are interested solely in the behaviour of the large 
particle, and hence the relevant quantity is the reduced density matrix, obtained by tracing 
over the environmental coordinates, 



«(.,.') ^n/''«..«(-.«".-'.^") 



(2.2) 



It is convenient to introduce the (non-unitary) reduced density matrix propagator, J, 
defined by the relation. 



Pt{x.,y) = / dxQdyo J{x,y,t\xo,yo,0) po{xo,yo) 



(2.3) 

Under the assumption that the initial density operator for the total system factorizes, 
P = P system ® Pbath^ ^^^ reduced density operator propagator is given by the path integral 
expression. 



i i i 

J{xf,yf,t\xQ,yQ,0) = I VxVy exp I -S[x] - -S[y] + -W[x,y] 



where 



h 



-Mx^ - V{x) 



h 



h 



S[x\ = dt 
and VF[a:(t), y{t)] is the Feynman- Vernon influence functional phase 

W[x{t),y{t)] = 



(2.4) 



(2.5) 



t rs 

ds / ds'[x{s) — y{s)] ri{s — s') [x{s') + y{s')] 
Jo 

+ i ds ds'[x{s)-y{s)]u{s-s') [x{s')-y{s')] (2.6) 

Jo Jo 



The kernels r]{s) and u{s) are defined by 

"°^ du 



Hs) 







TV 



liuj) coth I — — - I cos ujs 



d 



vis) =-7(.) 



(2.7) 
(2.8) 



Here 

r°^duj I(uj) , , 

7(s) = / L^cosc^s (2.9) 



and /(a;) is the spectral density 



n 



ImnOJn 



I{UJ) =\^b{uJ-UJn) ^_ " (2.10) 



The kernel z/(s) contributes a phase to the path integral (2.4), effectively modifying the 
action of the distinguished system. It leads to dissipation and frequency renormalization 
in the effective equations of motion. The kernel r\(^s) damps contributions from differing 
values of x and y. It is responsible for noise (and also for the process of decoherence 
discussed elsewhere). These two kernels are completely determined once a form for the 
spectral density (2.10) has been specified. A convenient class of choices is, 

I{uj) = M^uj (1)'"^ exp (^-^^ (2.11) 

Here, A is a cut-off, which will generally be taken to be very large, sometimes infinite, and 
a) is a frequency scale, which may be taken to be A. We will concentrate almost entirely on 
the case of the ohmic environment, s = 1, with occassional reference to the supra-ohmic 
and subohmic cases, s > 1 and s < 1 respectively. 

In the ohmic case, (2.8) is, 

7(s) = M7^ exp (-]A-'^sA (2.12) 

27r2 V 4 y 

and thus when A is very large, 

7(s) ^ M-f5{s) (2.13) 

We will work in this limit, unless otherwise stated. It should be noted, however, that 
this limit sometimes leads to a violation of positivity of the density operator at very short 
timescales (of order A^ ) [21]. We will say more about this later. 

In the limit in which (2.13) holds, the real part of the influence functional phase (2.6) 
may be written, 

ReW = ds M'y{x — y){x — y) 

Jo 

t 

ds M7 (5(0) {x^ - y^) + M7(x(0)2 - y{of) (2.14) 



6 



Here, the 5(0) is understood in terms of (2.12) at s = for large A. It is easily seen that 
the terms involving it may be absorbed by defining a renormalized potential in the path 
integral (2.4), 

Vjiix) = V{x) - M'yS{0)x'^ (2.15) 

For a harmonic oscillator of frequency ujq, considered below, we therefore define a renor- 
malized frequency, ujji, with uj^ = cUq — 27^(0). The terms involving the end-points a;(0), 
2/(0) are clearly negligible and will be dropped. 

The noise kernel (2.7) is non-local for large A, except in the so-called Fokker-PLanck 
limit, kT » HA, in which case one has 

.(>) ^ ^ m (2.16) 

An evolution equation for p may be derived. Its expected most general form is, 
- mrmx - J,) (^1^ - 1^ j - iT(t)h(t)(x - yfp 

+ ftrw/a)(x-.)(| + |) (2.17) 

The coefficients r(t), /(t), h{t), are in general rather complicated non-local functions of 
time, and appear to be known only in two particular cases. Firstly, explicit expressions 
for them may be found in Ref. [20] for the case of linear systems. Secondly, in the Fokker- 
Planck limit, for any potential V{x), one has [4], 

9MkT 
r(t)=7, h{t) = ^—, f{t)=0 (2.18) 

The propagator J may be evaluated exactly for the case of the simple harmonic oscil- 
lator, V{x) = ^MuqX . We will also be interested in the free particle, uq = 0. Introducing 
X = x + y,^ = x — y,it may be shown that [4,20], 

J{Xj, Cf, t\XQ, eo, 0) = — exp i-S - I J (2.19) 

where 

s = k{t)Xjij + k{t)XQiQ - L{t)XQij - N{t)XfiQ (2.20) 

and 

(/> = A{t)ef + B{t)ifiQ + C{t)il (2.21) 

Explicit expressions for the coefficients A, S, C and K, K, L, A^ are given in the appendix. 
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3. THE WIGNER FUNCTION PROPAGATOR 

Instead of the density operator, it is often convenient to work with the Wigner function, 
defined by 

W{p,q) = l^j di e-iP^ p(g + l^,g- 1^) (3.1) 

Its inverse is, 

p{x, y) = j dp eTA-~v) W{p, ^) (3.2) 

(See Refs. [22,23] for properties of the Wigner function.) 

The evolution of the reduced density operator is described by the (generally non- 
unitary) propagator, J, Eq.(2.4). Correspondingly, one may introduced the Wigner func- 
tion propagator, K{p, q, t\pQ, qq, 0), defined by 

K{p, q, t\po, go, 0) = ^ / ^^^^0 e-^P^J{q+^^, g-^^, t|go+^^0, 90-^^0, 0) e^P"^^> (3.3) 



and one has, 



Wip, Q.t) = J dpodqo Kip, q, t\po, go, 0) W{p, g, 0) (3.4) 



By performing the Wigner transform of the evolution equation for p, Eq.(2.17), one 
may derive an analogous equation for the Wigner function. It is 

dw p dw ,,, , ,dw ^^, , a , ,,,, ^^, ,,,,aV ,^, ,,, ,aV 

It follows that the Wigner function propagator also obeys this equation, together with the 
initial conditions, 

K{p, g, 0|po, ?0, 0) = 5{p - Po) S{q - go) (3.6) 

Consider the right-hand side of the Wigner function evolution equation (3.5). The first 
two terms are a Liouville operator, and if the equation had only these terms, W would 
evolve along the classical flow in phase space. The third term is the dissipative term. 
It modifies the flow along which W evolves, and also causes a contraction of each area 
element. The fourth and fifth terms are diffusive terms, and produce an expansion of each 



area element. They are responsible for noise/fluctuations, and also for the destruction 
of interference, by erasing the structure of the Wigner function on small scales. The 
last term, the power series, together with the first two terms make up the unitary part 
of the evolution. Hence, up to corrections of order h\ unitary evolution corresponds 
to approximately classical evolution of the Wigner function. It is partly for this reason 
that the evolution of a quantum system is most conveniently undertaken in the Wigner 
representation. The higher corrections can often be argued to be negligible, e.g., if the 
Wigner function does not develop too much detailed structure on small scales. There are, 
however, important examples where they cannot be neglected, e.g., in chaotic systems [24]. 

To illustrate the simplicity of evolution in the Wigner representation, we write down 
the solution for the propagator, in the case in which the diffusive terms, dissipative terms 
and higher order terms in h are ignored. It is, 

K{p, q, t\po, qo, 0) = 5{p - /^) 5{q - Q^') (3.7) 

where p = p {po, qo, t) and q = q {po, qo, t) are the solutions to the classical equations 
of motion, 

P=-V'(q) (3.8) 

satisfying the initial conditions p = pQ, q = qQ at t = 0. 

For the case of linear systems coupled to an environment, described in the previous 
section, the density operator propagator is given explicitly by (2.19)-(2.21). The Wigner 
function propagator may be computed, yielding, 

K{p, q, t\po, qo, 0) = -^ 

7ih{AAC-B^)2 

X e^p^-aip-p'^^f-Piq-q^V-^P-p'^-Q"^)) (3-10) 

where 

" = H4AC- B^) ^'-''^ 

^ _ AjAN^ + BNK + CK^) 

^~ H4AC - S2) ^^-^^^ 



2{NB + 2CK) 
' h{4AC - S2) 
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(3.13) 



and g , p = Mcf^^ are the solution to the classical field equations with dissipation and 
with a renormalized frequency, 

q + 2-iq + uj\q = Q (3.14) 

matching pQ, qQ oX t = 0. Explicitly, 

,^'^1^4,0 (3.15) 

p'^ = fpo-2(LiV-KK)| (3.16) 

Now for what follows, it is important that the transformation from po? 90 ^o p , Q' , 

defined by Eqs.(3.15), (3,16) is not a canonical transformation, since the Jacobian of the 

transformation is 

Q(^pd^qcl-^ _ L 

d{po,qo) N ^ • ^ 

This is because the classical evolution is dissipative, which tends to cause phase space cells 
to shrink. This shrinking effect can be compensated for by a scaling transformation of 
p", q^ . In particular, the variables, 

p' = \-^p^\ q' = X^^q^l (3.18) 

where A = {L/N)^, are related to pQ, q^ by a canonical transformation. Using these new 
variables, the propagation of the Wigner function may be written, 

W{p, q,t)= I dpodqo 



n(4AC-S2)2 
X exp (-a(p - Xp'f - (3{q - Xq'f - e(p - \p'){q - \q')) W{pQ, qo, 0) (3.19) 

Now one may perform the canonical transformation of integration variables from pQ, qq ^o 
p'q', 



W{p, q,t) = / dp dq 



n{4AC-B^)2 
X exp f-a{p - Xp'f - I3{q - Xq'f - e{p - Xp'){q - Xq')) W'{p', q', 0) (3.20) 



where 



W^'(p',Q',0) = Ti^(po,?0,0) (3.21) 
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Now the point is that the quantity W defined by this transformation is still a Wigner 
function, i.e., it is the Wigner transform of density operator, p', say. In fact, it is readily 
shown that p' is related to the original p by a unitary transformation. This would not be 
true of the transformation from pQ, qq to p , q . 

Eq.(3.20) is the main result of this section: a simple expression for the evolution of the 
Wigner function for linear systems from an arbitrary initial state. The (non-dissipative 
part of the) classical evolution has been absorbed into a canonical transformation of the 
initial state, and the effects of dissipation and diffusion are contained in the coefficients, 

a, (3, e, A. 

It is sometimes convenient to write the Wigner function propagator in the alternative 
form, 

K{p, q, t\po, qo, 0) = -^ 

7rh{4AC-B^)2 

X exp (-p(po - Pof - ^{qo - Qof - (j{po - Po){qo - Qq)) (3.22) 

where 

^ = H4AC- B^) ^'-''^ 



A{AK^ + ELK + CL 

h{4AC - S2) 
2{LB + 2Ak) 



(3.24) 

^ " n{4AC - S2) (^-^^^ 

Pq , Qq are again solutions to the classical field equations with dissipation, but now match- 
ing the final conditions p, q at time t: 

«S' = -^ + f. (3.26) 

p'o' = jP-2(LN-Kk)^ (3,27) 

We will use these expressions in the following sections. 



4. GENERALIZED UNCERTAINTY RELATIONS 

We now show how the results of the previous section may be used to derive generalized 
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uncertainty relations for quantum Brownian motion models, for linear systems. 

4(A). A Lower Bound on the Uncertainty Function 

From the Wigner function propagator, one may obtain expressions for the distributions 
of p and q at time t: 

piQ.Q^i) = dp W{p,q,t) 



dpodqo^-T exp I -il_^ 1 W'{p\ q' , 0) (4.1) 

(27rs2)2 V 2"? 



p{p^P^t) = dq W{p,q,t) 

1 

(27rs^)2 \^ ^Sp 



f dpodqo^-T exp (- (P y ) w'{p\ q\ 0) (4.2) 



where 



2 _ 2a 2 _ 2/? 

'^- (4a(3-e^y 'P ~ (4«/3 - e2) ^^'^^ 



From these results, it is straightforward to compute the variances of p and q at any time t: 

(Aqtf = X^Aq'f + si (4.4) 

(Apt)2 = A2(Ap')V«J (4.5) 

where (Ag') , (Ap') denote the variances of q and p in the canonically transformed initial 
Wigner function (3.21), or equivalently, in the unitarily transformed initial density operator 
Pq. (The variances are generally not invariant under such transformations of the state.) 

We have the usual uncertainty principle, 

{Aq'){Ap') > ^ (4.6) 

with equality if and only if p' is a coherent state, i.e., if and only if 

VVV.,'.0)^^e.p(-Hif:)!/-5^| ,4.7) 
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(where Aq is, so far, an arbitrary parameter). It follows that the uncertainty function at 
time t satisfies, 

Ut ^ {AqtfiAptf > ^ + X{Aq'fsl + ^^ + slsl (4.8) 

with equality if and only if the initial state is given by (4.7). Since Aq' is arbitrary, we 
may minimize over it. The minimum is at 

(Aq'f = ^ (4.9) 

and inserting the minimum value in (4.7), we obtain 

iAqt){Apt) > — + SpSq (4.10) 

In terms of the coefficents of the non-unitary propagator, (3.11)-(3.13), this reads, 

{Aqt){Apt) >^ + ^{ACN^ + BCNK + C^K^f (4.11) 

This the first main result of this section: the uncertainty at time t satisfies the generalized 
uncertainty relation (4.11), with equality if and only if the initial state is given by (4.7) 
with Aq' given by (4.9). 

Consider now the conditions for equality. From (4.7) and (3.21), the minimizing initial 
Wigner function is 



W{po, go, 0) = ;^ exp \--^—{Kpo - 2{LN - KK)qQY - —^^{-pQ + Kqo) 

(4.12) 
Inverting the Wigner transform, one finds that the minimizing state is a pure state, 

^(x)^(—^ exp(-{a + ib)x^) (4.13) 

This is a so-called correlated coherent state [18], discussed below. The coefficients a and b 
may be computed from (4.12) but will not be needed here so are not given. 

It is important to note that the initial state (4.13) minimizing the uncertainty at 
time t is different for each moment of time. That is, the lower bound (4.11), is not the 
time evolution of a particular initial state. It is actually an envelope. The initial state 
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(4.13), which depends on t, will minimize the uncertainty at time t but generally not at 
other times. This is because, although the effect of diffusion is generally to increase the 
uncertainty as times goes on, there are competing effects that may reduce it. In particular, 
wave packet reassembly (the time reverse of wave packet spreading) and dissipation may 
cause the uncertainty to decrease. This point is discussed further in Ref.[12], where an 
analagous lower bound on the Wehrl entropy (discussed below) was proved. 

4(B). A Sharper Lower Bound 

As stated, the lower bound on the uncertainty (4.11) is sharp, in the sense that at each 
moment of time there exists an initial state at which the lower bound is achieved. But it 
is only an envelope because the the initial state achieving it is different for each moment 
of time. It turns out that an improvement of this result is possible. 

First, we note that the the conventional uncertainty relation, ApAq > h/2, may be 
generalized to the stronger result [18], 

A^^{Apf{Aqf-C^pq>^ (4.14) 

where Cpq is the correlation between p and q, 

1 

2 

dpdq pq W{p, q) — {p){q) (4-15) 



Cpq =- {ApAq + AqAp) 



where Aq = q— (q), Ap — p = (p). The quantity A defined by (4.14) is essentially the area 
enclosed by the 1 — a contour of the Wigner function. Equality in (4.14) is achieved when 
the state is a correlated coherent state, 

„2 



*(^) = .. O.MA exp I -—2 



ir 
1 



(1 — r2)2 



+ ^-l(a2 + |«|2), (4.16) 



where ?7,r are real parameters (r = Cpq/ {ApAq) and r] = Aq), and a is a complex 
parameter. 

We have, however, already seen a correlated coherent state - it is the state minimizing 
the uncertainty relation (4.11). It is therefore plausible that a stronger version of (4.11) 
might be possible in terms of the Wigner function area (4.14). Indeed, computing the 
correlation from (3.19), one finds the remarkably simple result. 



L2_ 2 fl^ 



i2 _ ^ ,2 , r^ /, A^ d2x 



^t =^^0 + Jp^i^AC-B^) (4.17) 
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This gives the Wigner function area A at an arbitrary time t, At^ in terms of its value at 
time t = 0, ^0- It is remarkable because A at time t depends on only its value at t = 0, 
and not on any of the other moments of the initial state. Now at t = 0, Aq satisfies the 
uncertainty principle (4.14). We thus obtain, for the Wigner function area at time t, the 
inequality 

Equality is achieved when the initial state belongs to the four-parameter family of corre- 
lated coherent states. Most importantly, the inequality is sharp at every moment of time 
for one and the same initial state. It is not an envelope. Eq.(4.18) is the second main 
result of this section. 

The simplicity of the time evolution in terms of the the Wigner area is easy to un- 
derstand. For linear systems, the area is preserved under unitary time evolution. For the 
non-unitary case considered here, therefore, (4.17) expressess entirely the effects of the en- 
vironment, with the effects of unitary evolution completely factored out. In this respect it 
is similar to the von Neumann entropy discussed in the next section. Indeed, for Gaussian 
states, the von Neumann entropy is a function only of the Wigner area. 

Some rather different, but not unrelated generalizations of the uncertainty principle 
may be found in Refs. [25,26]. 



5. EXPLICIT FORMS OF THE LOWER BOUNDS 

We now give the explicit forms of the lower bounds (4.11) and (4.18), in a variety of 
situtations. There are a very large number of different situations and regimes that our 
results may cover, depending on the choice of t, T, 7 and uj (harmonic oscillator in the 
under- and over-damped case, inverted harmonic oscillator, free particle). We will not 
present here a completely exhaustive search of the parameter space, but concentrate on 
the cases that contain some interesting physical results. In particular, most of the following 
results are the underdamped case. The cases not covered here are readily derived, should 
the reader be interested. 

Consider first the case of short times, t << 7^ . Then for high temperatures, (4.11) 
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and (4.18) become, respectively, 

These relations are in fact valid for any potential. They represent the initial growth of 
fluctuations, starting from the purely quantum fluctuations at t = 0. 

The 7t term in each expression indicates an initial decrease of fluctuations, in apparent 
violation of the uncertainty principle. This violation occurs, in Eq.(5.1) for example, on a 
timescale less than 7i/kT. This is because these expressions have been derived by taking 
the inflnite cut-off limit in Eq.(2.13), which as previously stated, can lead to violations 
of positivity of the density operator on time scales less than A^ [21]. The expressions 
are therefore not valid for t < A^ . The high temperature limit used in deriving (5.1) 
means kT » hA. Combining these inequalities, we see that (5.1) is valid only for times 
t > A~ >> h/kTj for which violations of the uncertainty principle will not arise. 

To illustrate explicitly that it is the inflnite cutoff limit that is responsible for violations 
of the uncertainty principle, we compute the lower bound (4.11) leaving the cutoff flnite in 
the expressions for the various coefficients in the propagator (2.19)-(2.21). The important 
point is to keep A flnite in Eq.(2.12). Expressions (A.1)-(A.4) and (A. 14) are no longer 
valid, and the correct expressions for these coefficients, in the short time limit, are given 
by Eqs. (A.28)-(A.32). One thus obtains, in the short time limit, the lower bound 

(Ag,)(Ap,) > | + ^t3 + ... (5.3) 

for which there is clearly no violation of the uncertainty principle. 

Ignoring the positivity-violating terms in (5.1), (5.2), the remaining terms give an indi- 
cation of the time scale on which the thermal fluctuations become important in comparison 
to the quantum ones. It is. 

As noted in Ref.[27], this is the timescale on which quantum fluctuations become compa- 
rable to Nyquist noise. 
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The timescale (5.5) differs from that derived in Refs. [10,12], where it was argued 
that thermal fluctuations become comparable to quantum ones on the timescale t^i ~ 
h /{M'fkTi ) (the timescale characteristic of decoherence [28]). Here £ is a length scale, 
arising from the particular way in which the fluctuations were measured in Refs. [10,12] 
{e.g., in Ref.[12] is was the width of a phase space projector). There is no length scale 
present in the method use to measure fluctuations here, hence a different result is obtained. 
The exact connection between these two results is not very clear. We note, however, that 
the time scale (5.4) will typically be much longer than the decoherence time scale. 

In the weak coupling regime, for all temperatures and times (except very short times), 
Eq.(4.18) becomes, using (A.18)-(A.20), 



2, ^2 r 



At> 



4 



}-'"^^^'^i^)^'-^-'''r 



(5.5) 



An expression may also be derived for arbitrary coupling and high temperatures, using 
(A.11)-(A.13), but this is rather complicated and will not be given. 

In the long-time limit, for any coupling and temperature. 





Ut = 


A> 


-I 




In the high temperature limit this becomes. 








Ut 


= A't 


k^T^ 
^R 


and 


in the weak coupling limit. 










Ut 


= A 


> 

- 4 





(5.6) 



(5.7) 



(5.8) 



6. INEQUALITES RELATING UNCERTAINTY 
TO ENTROPY AND LINEAR ENTROPY 

The von Neumann entropy of the reduced density operator, 

S[p] = -Tr{plnp) (6.1) 
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is often discussed in the present context, in association with uncertainty, decoherence, and 
correlations of the distinguished systems with its environment [14,15,17,29,13]. Zurek, Paz, 
and Habib for example, looked for classes of initial states, which under evolution according 
to a master equation of the form (2.17), generate the least amount of entropy at time 
t. They regarded such states as the most stable under evolution in the presence of an 
environment. They argued that the initial states doing the job are coherent states, at least 
approximately. 

One of the reasons for looking at the von Neumann entropy is that it is constant for 
unitary evolution, thus for open systems such as those considered here, it is principally 
a measure of environmentally induced effects. The Wigner area considered in Section IV 
also has this property. The results of Section IV thus agree with the claims of Zurek et al. 
(except that it is really the correlated coherent states, rather than the ordinary coherent 
states which are the most stable, although this distinction is not very important). It would 
be useful to find a connection between these two measures of stability or fluctuations. 

6(A). Uncertainty vs. von Neumann Entropy 

The connection between uncertainty and the von Neumann entropy may be found 
indirectly, by considering flrst of all the phase space distribution. 



^{p,q) = {z\p\z) (6.2) 



where 



{x\z) = (x|p, ,) = I J-2 j exp [-^-^j^ + ^P^j (6.3) 

are the standard coherent states. The function n{p, q) is normalized according to 

/^.(.P)^l. (6.4) 

It is readily shown that /u(p, q) is also equal to 



2al 2a2 / '' P' 



nip, q) =2 I dpdq exp ( ^^^ ^^-^ — ) Wp{p , q ) (6.5) 
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where Wp{p, q) is the Wigner function of p and apaq = ^%. Eq.(6.5) is sometimes known 
as the Husimi distribution [30], and is positive even though the Wigner function is not in 
general (see also Ref.[31]). 

There exists an information-theoretic measure of the uncertainty or spread in phase 
space contained in the distribution (6.5), namely the so-called Wehrl entropy [32] 

/(P, Q) = - I ^ ^(p, q) ln^(p, q) (6.6) 

This is the Shannon information of (6.2). /(P, Q) is large for spread out distributions, and 
small for very concentrated ones. Because of the uncertainty principle, one would expect 
a limit on the degree to which p{p, q) may be concentrated about a small region of phase 
space, and hence a lower bound on (6.6). Indeed, a non-trivial theorem due to Lieb shows 
that, 

/(p,g) > 1 (6.7) 

with equality if and only if p is the density matrix of a coherent state, \z'){z'\ (Ref.[33]). 
In Ref.[12], the Wehrl entropy was used as a measure of both quantum and thermal fluc- 
tuations, and a lower bound analagous to (4.11) was derived. 

The reason for studying this quantity is that it provides the link between the von 
Neumann entropy and the uncertainty measures U and A. On the one hand, an elementary 
property of Shannon information is 

/(P,Q) < ln(|(detK)5) (6.8) 

where K is the 2x2 covariance matrix of the distribution p{p, q) [34]. Equality holds if 
and only if p{j>, q) is a Gaussian. From (6.5), one has 

detK = ((Aq)2 + af) [{Apf + a2) - C^^ (6.9) 

where Aq, Ap are the quantum-mechanical variances of the state p. On the other hand, 
an elementary property of the Wehrl entropy, following from the concavity of Shannon 
information is [32], 

/(p,g) > s[p] (6.10) 

Hence combining the upper and lower limits (6.8), (6.10), one obtains, 

{{Aqf + a2) ((Ap)2 + a^) - C^^ > h^e^^^^-^^ (6.11) 
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Finally, since the parameter aq is arbitrary (and apaq = ^/2), we may minimize the 
left-hand side over it, with the result, 

(ApAg+Jn)2-C2^ > nV(^-l) (6.12) 

This is the exact form of the connection between the uncertainty and entropy for a general 
mixed state p. 

In the regime where quantum fluctations are more significant than thermal ones, it is 
appropriate to use the lower bound (6.7) rather than (6.10) (since S[p\ goes to zero if the 
state is pure), and this is formally achieved by setting 5" = 1 in (6.12). One then deduces 
the usual uncertainty principle from (6.12) (although not the generalized version including 
the correlation, (4.14)). 

In the regime where thermal (or environmentally-induced) fluctuations are dominant, 

one would expect ApAg >> ^/2 and 5" >> 1, and (6.12) then gives 

ApAo A q , ^ 

— V- ^ > ^ > e^ 6.13 

Ti n 

This is the simplest form of the connection between the uncertainty and the von Neumann 
entropy: the entropy is bounded from above by the logarithm of the number of phase space 
cells the state occupies. 

To see how sharp these equalities can be, consider the case of a Gaussian Wigner 
function. It may be shown that the von Neumann entropy of a Gaussian is 

where n = 2A/h and A is the area of the Gaussian Wigner function (see Ref.[13] for 
example). For large A, 

S[p] ^ In (^) (6.15) 

and hence we have equality in (6.13). 

6(B). von Neumann Entropy vs. Linear Entropy 

In practice, the discussions of Zurek et a/., and indeed many discussions of von Neu- 
mann entropy, often concern the so-called linear entropy, 

Sl = 1- Trp2 (6.16) 
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How is this related to von Neumann entropy? Let the density operator be 

p = ^ Pn \n){n\ (6.17) 

n 

We use Jensen's inequahty (see, Ref.[34], for example), which states that if / is a convex 
function, and X a random variable, 

{f{X)) > f{{X)) (6.18) 

where (■ ■ ■) denotes the mean over X . Take / to be the exponential function, and X to be 
In pn- Then, 

^pI > expi^pnlnpnj (6.19) 

n \ n / 

and thus 

Trp2 > e~^^P^ (6.20) 



In terms of the linear entropy. 



Sl = 1- Trp2 > 1 - e"^ (6.21) 



Equality in (6.21) is reached for pure states, when S = Si = 0, and is for very mixed 
states, when 5" is very large and 5'^ ~ 1. 



7. LONG-TIME LIMITS 

One of the particularly interesting questions for non-equilibrium systems of the type 
considered here is whether they settle down to a unique state after a long period of time. 
It turns out to be particularly straightforward to answer this question for linear systems, 
using the Wigner function propagator described in Section III. To see this, note that 
combining (3.4) and (3.10), we have the expression, 

W{p, q, t) = -^ / dpodqo W{po, qq, 0) 

nh{4AC-B^)2 J 

X exp (-a(p - /^)2 - (3{q - q^f - e{p - /^) {q - q''^)) (7.1) 

21 



where, recall, p^\q^^ are the value of the classical solutions (with dissipation) at time t 
matching the initial datapo? Qo Sitt = 0. This expression allows us to compute the moments 
at any time t in terms of the initial moments. By computing the long-time limits of these 
moments, the form of the long-time limit of the Wigner function may be obtained, since 
it is completely determined by its moments. 

For the harmonic oscillator, simplifications occur. Consider first the case of the har- 
monic oscillator in an ohmic envrionment. One has p = Mq and q satisfies Eq.(3.14). 
Either from Eq.(3.14), or from the explicit solution (3.15), (3.16), it is easily seen that 
p — > and q ^ in the long time limit. All dependence on pQ and go drops out of the 
exponential in (7.1), and one obtains the following expression for the the asymptotic value 
of the Wigner function, 

W^{p,q) = -^exp -ap - (3q - epq) (7.2) 

7rh{4AC-B^)2 ^ ^ 

The coefficients a, /3, e are given by the long-time limits of (3.11)-(3.13). Using results 

(A.21)-(A.23) in the appendix, one finds 






e =0 (7.4) 

with Fi and F2 given by (A. 24), (A. 25). These relations represent the exact form of the 
long-time limit in the ohmic case. 

It is useful to compare the result (7.2) with the Wigner function of the harmonic 
oscillator in a thermal state: 

(7.5) 



»Vfe .) = ;^ tanh (II) exp (- tanh (5^) 



Mujp 9 1 9 



2kT MhujR 

The long-time limit (7.2) coincides with the thermal Wigner function (7.5) in the high- 
temperature limit, and in the weak coupling limit, as may be seen from Eqs.(A.26), (A. 27) 
in the appendix. Note, however, that the asymptotic state is not always a thermal state. 

It is also of interest to compute the uncertainty in q in the long time limit. It is, 

- ^ io '" "'^ [^) M (,|_, 2)2 + 4^2,2 (7-6) 
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This is in agreement with the fluctuation-dissipation relation [35]. Eq.(7.6) was given in 
Ref. [4] , for the special case of a Gaussian initial state, but not surprisingly coincides with 
this more general result because every initial state goes to a Gaussian. 

For the harmonic oscillator with a non-ohmic environment, it may be shown, using the 
more general treatment of Ref. [20], that q in the Wigner function propagator in (7.1) is 
the solution to the integro-differential equation, 

rt 
q{t) + ulq{t) - / ds 7]{t - s) q{s) = (7.7) 

Jo 

where r](t — s) is dissipation kernel (2.10) for an arbitrary spectral density I{uj). We 
conjecture that there exists a large class of spectral densities for which all solutions to 
(7.7) will satisfy q -^ 0, p ^ in the long-time limit, although we have not been able 
to prove this. If this conjecture is true, then the long-time limit of the Wigner function 
will again be of the form (7.2) (although the coefficients ct, /3, e will not necessarily be the 
same) . 

These results, including the above conjecture, are consistent with the proof in Ref. [11], 
that the Wigner function of a single member of a chain of coupled harmonic oscillators 
tends towards a Gaussian Wigner function in the long-time limit, under certain reasonable 
conditions on the environment. 



8. GENERAL POTENTIALS AND THE FREE PARTICLE 

All of the results we have described so far concern the specific case of linear systems. 
It is clearly of interest to generalize our considerations to other types of potential, and this 
we now do. 

We are interested in the propagator (2.4) for arbitrary potentials V{x). It is generally 
not possible to evaluate this propagator for arbitrary V{x), so some approximations need 
to be made. Again it is most convenient to work with the Wigner function evolution 
equation (3.5). The coefficients of the dissipative and diffusive terms in (3.5) are known 
only in either the linear case treated earlier, or in the high-temperature limit, so we consider 
the latter. We will also neglect the infinite power series in derivatives of W and V, and 
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hence the equation becomes, 

We expect that the neglect of higher derivative terms in W and V is vahd when the initial 
Wigner function does not possess, or develop, too much detailed structure in regions of 
size h or less. It is argued in Ref.[24] that the diffusive terms rapidly smooth the Wigner 
function on small scales, and thus (8.1) is plausibly valid at long times for any initial state. 

As in the linear case, it is possible to associate a Wigner function propagator (3.4), 
with (8.1), although it is not possible to calculate it explicitly. We have explored the 
possibility that the Wigner function propagator is given approximately by an expression 
of the form (3.10), where p^/ q,d ^^^ ^j^g solutions to the dissipative equations of motion 
with potential V{x), and the coefficients a, /?, e are determined by substituting (3.10) 
into (8.1). Unfortunately we have not been able to justify this approximation, and as a 
consequence, we have not extended our results on the time-dependent lower bounds on 
uncertainty functions derived in Section IV to arbitrary potentials. We will therefore work 
directly with Eq.(8.1), and consider only the question of the long-time limit. 

8(A). Long Time Limits for General Potentials 

Eq.(8.1) is the Kramers equation. Its properties are described in some detail in the 
book by Risken [36] , and we shall make use of those results here. Many of the results in 
Risken concern non-negative phase space distribution functions, but we have been careful 
not to assume that here, because the Wigner function is of course not always positive. 

It is easy to see that the master equation possesses the stationary solution, 

.2 



H'fe«)-i*e.p(-^-l|M) (8.2) 



where A^ is a normalization factor. This will be an admissable solution, ie., a Wigner 
function, only if the potential is such that exp (—Vji{q)/kT) is normalizable. This requires 
^i?('?) ^ oo as g ^ ±oo faster than In Igj. In that case, the stationary distribution is then 
the Wigner transform of a thermal state, p = Z e^ ' where Z = Tr(e^ ' ), for 
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large temperaturesT . 



The master equation may be written 

dW 



CW (8.3) 



dt 

where C is an operator deducible from (8.1), and acts on the set of square integrable 
functions on phase space, with the inner product between functions defined by 

(/, 9)= j dpdq f{p, q) g{p, q) (8.4) 

The general solution to the master equation may then be written 

W{p, q,t) = Y, Cn e-^-* Mp, q) (8-5) 

n 

where the c^ are constants, and (/)n are the eigenfunctions of £, with eigenvalue A^^, 

C(f)n = \n4>n (8.6) 

It may be shown that if the master equation admits a stationary solution, i.e., if one of the 
eigenvalues is zero, then the real part of the remaining eigenvalues is strictly positive [36]. 
It immediately follows that the general solution will tend towards the unique stationary 
solution (8.2) in the limit t -^ 00. Hence every initial state tends to the thermal state in 
the long time limit. 

For potentials which are not bounded from below, the long-time limit will often depend 
quite sensitively on the initial state. Indeed, the inverted oscillator was recently used by 
Zurek and Paz as a prototype model for chaotic systems [24]. 

8(B). Long Time Limit for the Damped Free Particle 

Finally, let us consider the case of the free particle, ie., Vji{q) = in Eq.(8.1). In 
this case, there is no stationary solution, because the operator C does not have any zero 
eigenvalues. To see this, let us look for eigenfunctions of the form 

Unk{P.q)=e'^'^fn{p) (8.7) 



T Actually, the only Wigner function that can be positive is a Gaussian one [23]. However, 
there is no contradiction with (8.2) because we are working with a distribution function 

that satisfies the full Wigner function equation only up to order fi , so the distribution 
function is not exactly a Wigner function 
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In Ref.[36], it is shown that the eigenvalue equation for a zero eigenvalue, Cu^}^ = 0, is 
satisfied only if /c = 0. Hence although a formal solution to the eigenvalue equation exists, 
it does not belong to the spectrum of C because it is independent of q and so will not be 
normalizable. 

Despite the absence of a stationary solution, it is still of interest to ask whether the 
one can say anything at all about the long time limit of the solution. Consider Eq.(7.1) 
for the damped free particle. For large t, one has 

1 . M7 1 



2MkT 2kTt 2kTt ' ' 

A = B = ^, c = HM^ (8.9) 



Introducing q = Qq + pQ/{2M'y), Eq.(7.1) may be written. 



where 



W{p, q, t) = -^ exp - a - — p 

Tvh{4AC-B^)2 V V 4/3; 

X I dqgoiq) exp\-/3(q + ^p-q] | (8.11) 

90iq) = J dpo W{po, q - ^, 0) (8.12) 

The integrand of (8.12) is still a Wigner function, since the shift in the q argument can be 
compensated for by a unitary transformation of the density operator. 

If we integrate out q, and noting that e /4/3 -^ for large t, we obtain 

dq W{p, q, t) = i r exp ( -^r^] (8.13) 



{2'KMkT)2 \ 2MkT^ 

Hence the distribution of momenta approaches a Boltzmann distribution for all initial 
states. The remaining question is, what we can say about the q distribution in the long 
time limit? 

Let y = q+jgP = g— 2M-y - Then the integral in the expression for the Wigner function 



^.11) is, 



g{y,t) = dq goiq) exp{-(3 iq + —p-qj j 



dq 90{q) exp ("^^^^ ) (8.14) 
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where D = kT/{2M'y). g{y,t) obeys the diffusion equation 

From (8.14) or (8.15), one may compute aU the moments of y at time t in terms of their 
initial moments. One has, for example, 

{y^)t = 12I^V ^ i2Dt{y\ + {y\ (8.16) 

{y% = 60I)V(y)o + 20Dt{y\ + {y\ (8.17) 

From the explicit expressions for the moments such as these, it is straightforward to show 
that the leading order behaviour of all moments of the form (y") as t ^ oo are correctly 
reproduced by the Gaussian distribution, 

»,,) ^ -1-, exp (-<^) (8,18) 

It follows that for large times, the Wigner function approaches the asymptotic form, 

W(p,,,t) ^ exp (-J^) exp (-m: (« - ^ - («)0 - m'] (8.19) 



2MkTj ^ \ 2kTt V 2M7 '^'" 2M-fJ 

Hence the distribution of p is a Boltzmann distribution, as noted above, and the q distri- 
bution, obtained by integrating out p, is peaked about the value q = {q)Q + (p)o/(2M7), 
which is the asymptotic value of q under classical evolution, starting from the initial values 

(?)o, {p)o- 

The result (8.19) can, however, be rather misleading, and has limited value as an 
approximation to the Wigner function at large times. To understand this, write the Wigner 
function for large times as 

W{p,q,t) = Wsip,q.t) + Wi{p,q,t) + ... (8.20) 

where Wg is the leading order approximation to the Wigner function for large t and Wi is 
the next to leading order term. Let us compare the free particle case considered here with 
the case in which there is a stationary solution. When a stationary solution exists, Wg 
is the stationary solution, and hence is independent of time. Furthermore, from Eq.(8.5), 
the next term Wi is proportional to e^ 1* where Ai > is the first eigenvalue. It follows 
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that all moments of (8.20) are given by their moments in Wg plus an exponential decaying 
term. Even if the moment of Wg vanishes, the correction given by Wi goes to zero for 
large t. All moments of (8.20) approach the stationary moments like e~ ^ . 

Now consider the free particle case. Here there is no stationary solution, so Wg is time 
dependent. The next correction Wi does not decay exponentially fast. There are certain 
moments (e.g., {{y — (y)) )) that are zero for Wg and take their leading contribution from 
Wi- Such moments grow with time, as may be seen from (8.16), (8.17), so unlike the 
stationary solution case, they are not well approximated by their moments in Wg. Also, 
even for the moments which are well-approximated by the moments in Wg for large t, the 
rate of approach to the regime in which that approximation is valid depends on the initial 
conditions. It does not proceed at a universal rate. 

Therefore, although Eq.(8.19) is the formal solution to the Wigner function equation 
for large t, it is not very useful. In practice it will be more useful to work directly with 
the equations for the moments. 



9. SUMMARY AND DISCUSSION 

We have studied the evolution of open quantum systems described by the evolution 
equation (2.17). We were concerned with two particular questions: generalized uncertainty 
relations for this class of non-equilibrium systems and long-time limits. Our results may 
be summarized as follows: 

1. For any linear system whose evolution is described by the propagator (2.19)-(2.21), 
the uncertainty U and the Wigner function area A have the sharp lower bounds (4.11), 
(4.18), respectively. These represent the least possible amount of noise the system must 
suffer after evolution for time t in the presence of an environment. These expressions 
are valid for all types of environment {i.e., for all choices of spectral density). Also 
worthy of note is the particularly simple expression (4.17) of the evolution of the Wigner 
function area A. 

2. For the particular case of the ohmic environment, the explicit form of the lower bounds 
is given in Section V. These explicit expressions give the comparative sizes of quantum 
and thermal fluctuations, generalizing the work of Hu and Zhang [10]. 
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3. For the linear systems considered here, these generahzed uncertainty relations achieve 
equality for Gaussian pure initial states of the form (4.16). Such states are therefore 
the ones that suffer the least amount of noise under evolution in the presence of an 
environment, substantiating the results of Zurek et al. [17,14,15]. 

4. The uncertainty is connected to the von Neumann entropy via the relation (6.12), and 
the entropy is connected to the linear entropy by the relation (6.21). 

5. For a harmonic oscillator in an ohmic environment, all initial states tend towards a 
Gaussian Wigner function in the long-time limit. It is the same as a thermal state in 
the high temperature limit or the weak coupling limit. For non-ohmic environments, 
the same result is plausibly true for a reasonable class of spectral densities, but this 
remains to be shown. 

6. For general potentials such that exp {—V^{q)/kT) is normalizable, the Wigner function 
tends towards a thermal state in the high temperature limit, for all initial states. For 
the free particle, the Wigner function tends to a Gaussian state, although this is not 
a very useful expression because it does not give a correct approximation to all the 
moments for large times. 

The reason we were able to prove the results (1) and (4) with such ease was our 
use of the Wigner function propagator, in terms of which the quantum evolution takes a 
particularly transparent form. We comment that the detailed methods used here could well 
be of use in related calculations. For example, it might be possible to discuss decoherence 
of arbitrary initial states using the Wigner function propagator derived in Section III. 
These and other related questions will be pursued in future publications. 
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APPENDIX 

In this appendix we give the exphcit forms of the coefficients A, B, C and K, K, L, N 
appearing in the exphcit expression for the propagator, (2.17). The foUowing results are 
taken from Caldeira and Leggett [4], and from Hu and Zhang [10], with minor elaborations 
and extensions. 



We first give the forms of the coefficients for the harmonic oscillator in the under- 

2 
R 

Then we have, 



damped case, w^ > 7. Let uj = cun — 7 . We work in the underdamped case, 7 < ujji. 



Also, 



where 



Similarly, 



where 



K{t) = --M-f + -Muj cot ujt, {A.l 

K{t) = +-M-f + -Moo cot ut, {A.2 

Lit) = , M.3 

^ ' 2smujt ' ^ 

Nit) = TT— — • (^-4 
Zsmujt 



A{t) = —^j^ djy expl-j^j u coth\--^j A^{t) (A5 

dr ds sin cur cos z/(r — s) sincus e''''^^^*'' {A.6 

Jo Jo 



e-27t ft 
AAt) 



sin^ ut 



M'y f°° ( v^\ ( Uv\ 

B{t) = ^y^ du exp (^--^j u coth i^^ B,{t) {A.7 



2e^^ f f ( ) 

Bi/{t) = — 7^ — I dr I dssirvuJT cos z/(r — s) sina;(t — s) e'*''-'^''''^-' (A. 8 

sin^^cut Jo Jo 

and 

C{t) = ^ j^dv exp (^-^^ z. coth (^^^ C,{t) (A.9 



where 

ft rt 

dr / 
sin^cut Jo Jo 
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Cj,(t) = — ^ — I dr I dss\nuj{t-T) cos z/(r - s) sincu(t - s) e^(^+*) (A.IO 

sin^cut Jo Jo 



We have included, for completeness, the explicit dependence on the cutoff in the expressions 
for A, B and C, and this is sometimes required, although we will generally work with the 
case A -^ oo. The integrals for Ay, By and Cy have been evaluated by Hu and Zhang [10], 
and we will make heavy use of their results. The remaining integrals over z/ to yield A, B 
and C cannot be carried out in general, but asymptotic expansions are possible in various 
regimes of interest, and these we now give. 

(Al). High- Temperature Limit 



In the much-studied high-temperature (Fokker- Planck) limit, one has coth(^z///cT) ^ 
kT/hu, and the integrals (A.5)-(A.10) may be evaluated exactly for any t, with the results. 



A = 
B = 



C 



MkT 
2% sin ujt 
MkTe'^ 
h sin ujt 

MkTe'^^^ 



1 



,-27* _ 



1 2 \ 7 

—^ ^(7 COS 2a;t + 0^7 sin 2a;t) ^ f 

7^ + 0;^ 7^ + UJ^ 

1 



-2^t 



-(1 - e"^^^) cosut + -^ ^ (7^(1 + e"2^^) cosc^t 



UJ 



2h sin ujt 
At short times, one has 

and for long times. 



_7^ -|- w^ 



- e-27* + 



7^ + cu^ 

-Fcu7(l-e"^'^*)sincc;t 
e-27i 



7^ -|- w^ 



(7 cos cut — a;7 sin 2a;t) 



A = B = C 



A 
B 
C 



MkT 
2^ sin cut 
MkTe^^ 
h siir' ut 
MkTe'^^^ 



1- 



7^ + cu^ 



2M-fkTt 

3h 



(7 COS 2c<;t + 0^7 sin 2c<;t) 



UJ 



7^ + cu^ 



COS cut + 



CU7 



7^ -|- cu^ 



sin cut 



cu 



2n sin^ cut (7^ + cu2) 



(All) 

(A.12) 
(A.13) 

(A.14) 

(A15) 
(A.16) 
(A.17) 



(A2). Weak Coupling Limit 



In the weak coupling regime, 7 << cu, but for arbitrary temperatures, one has, from 
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Hu and Zhang, 



. Muj ( huj 

A — K — cotn 

4 sin^ ut 

B = ^ — coth 



2 sin^ ut 

G = ;5 — cotn 

4 sin^ ut 



2kT 

huj 

2kf 
2kT 



1 - ^ sin 2ujt - e~2^* 



1 1 - 

— sin ut — cos ut + (— sin cut + cos cut) e" 

u u 



-2jt 



1- (l + -sin2cut)e 

u 



-2jt 



(A18) 
(A19) 
(A.20) 



where terms of order "f /u have been neglected. The long-time limits of these expressions 
are easily seen. 

(A3). Long-Time Limit for Arbitrary 
Temperature and Coupling 



It is also possible to determine the exact form of the long-time limits of A, B and C, 
without assuming high temperature or weak coupling. This is necessary in order to give 
a completely general statement about the long-time limits of arbitrary initial states, as in 
Section VII. From Hu and Zhang, they are 



A 



M 



2u'j^ sin^ ut 



2 2 T 

(u COS ut — '-yu sin 2c<;t) (Fi -\ — F2) 

u 



+u^sm'^ut {{1 + 



2Y 



u^ 



7 



u 



F2 



B 



C 



Me^^ 2 7 

(70; sin cut — cu coscut)(FiH F2) 

u 



cu D sin^ cut 

Mcu^e^T* 
2cu^ sm^ cut cu 



7 



where 



Ft = — / du u coth — — - 
27V Jo \2kT 

t'2 = — / au u coth 

27r Jo 



7 



+ 



7 



V2A;tJ 



7^ + (cu + z/)^ 7^ + (cu - u)' 
u + u u — u 



+ 



(A21) 

(A.22) 

(A.23) 

(A.24) 
(A.25) 



^72 + (cu + z/)2 -^2_^(^_^)2^ 

The integrals Fi and F2 do not appear to be exactly soluble, but can be evaluated in the 
high temperature and weak couling regimes. In the high temperature regime, 

kT 
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Fi 



F2^0 



(A.26) 



and the results (A.11)-(A.13) are recovered. In the weak couphng regime, 

Fl«i-oth(|^), F,«0 (A27) 
and the results (A.18)-(A.20) are recovered. 

(A4). Very Short Time Limit vi^ith Finite Cutoff 

For finite cutoff, and in the very time limit, t << A^ , and at high temperature, 
kT » huj, one has 

A = 2B = C = ^^^t' (A.28) 

and 

K = — H '- — r (A.29) 

2t 457r ^ ^ 

K = — + -^t^ (A.30) 

M IIM7A3 o , , , 

2t ISOtt ^ ' 

iV = ^ — r A.32 

2t ISOtt ^ ^ 
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